Spin and Valley dependent analysis of the two-dimensional low-density electron 

system in Si-MOSFETS. 
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The 2-D electron system (2DES) in Si metal-oxide field-effect transistors (MOSFETS) consists of 
two distinct electron fluids interacting with each other. We calculate the total energy as a function 
of the density n, and the spin polarization £ in the strongly-correlated low-density regime, using a 
classical mapping to a hypernetted-chain (CHNC) equation inclusive of bridge terms. Here the ten 
distribution functions, arising from spin and valley indices, are self-consistently calculated to obtain 
the total free energy, the chemical potential, the compressibility and the spin susceptibility. The 
T = results are compared with the 2-valley Quantum Monte Carlo (QMC) data of Conti et al. 
(at T = 0, f = 0) and found to be in excellent agreement. However, unlike in the one-valley 2DES, 
it is shown that the unpolarized phase is always the stable phase in the 2-valley system, right up 
to Wigner Crystallization at r 3 — 42. This leads to the insensitivity of g* to the spin polarization 
and to the density. The compressibility and the spin-susceptibility enhancement calculated from the 
free energy confirm the validity of a simple approach to the two-valley response based on coupled- 
mode formation. This enables the use of the usual (single- valley) exchange-correlation functionals in 
quantum calculations of MOSFET properties provided mode-coupling effects are taken into account. 
The enhancement of the spin susceptibilty calculated from the coupled-valley response and directly 
from the 2-valley energies is discussed. The three methods, QMC, CHNC, and Coupled-mode theory 
agree closely. Our results contain no ad hoc fit parameters. They agree with experiments and do 
not invoke impurity effects or metal-insulator transition phenomenology. 
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I. INTRODUCTION. 

The 2-D electron systems in GaAs-based structures as 
well as those found in Si MOSFETs are being intensely 
studied owing to the accessibility of a wide range of elec- 
tron densities under controlled conditions, leading to a 
wealth of the experimental observations The nature 
of the physics depends very much on the "coupling pa- 
rameter" r = (potential energy) /(kinetic energy). The T 
for the 2DES at T = 0, and the mean density n is equal 
to the mean-disk radius r s = (7m) -1 / 2 per electron, usu- 
ally expressed in effective atomic units which depend on 
the bandstructure mass and "background" dielectric con- 
stant. The 2DES in GaAs/AlAs structures will be called 
a simple 2DES or one- valley 2DES to distinguish it from 
the two- valley system found in Si-MOSFETS. The inver- 
sion layer adjacent to an oxide layer grown on the Si- 
(001) surface contains two equivalent valleys which host 
the two valleys of the electron fluid. Various aspects of 
such multi- valley systems were studied by Sham and 
Nakayama, Rasolt et al, and others, mainly in the high 
density limit. The simple 2DES is also a two component 
system because of the two spin species, while the 2-valley 
system involves 4-components and ten pair-distribution 
functions (PDFs). 

In a recent study of the effective mass m* and the effec- 
tive Lande-g factor of the 2DES, the Coulomb coupling 
between the electrons of the two valleys was shown to 
have a dramatic effect at low densities when the coupling 



becomes large^. In effect, the elementary excitations of 
the two fluids in the two valleys interact to form cou- 
pled modes, giving rise to new effects. It was shown 
experimentally!^] that m* g* rises rapidly with decreas- 
ing density in Si-MOSFETS, and that this rise is due to 
a dramatic increase in m* , independent of the spin polar- 
ization, while g* remains essentially constant. Calcula- 
tions for the Si system which account for the inter-valley 
Coulomb coupling quantitatively predict Q the sharp in- 
crease in m*g*. It was also shown that g* remained essen- 
tially constant, in strong contrast to the behaviour found 
for the simple one valley 2DES0. The effective mass 
was also shown to be practically independent of the spin 
polarization £, in excellent quantitative agreement with 
the data of of Shashkin et al.0 for Si-MOSFETS. This 
result leads to the view that the rapid rise in m* near 
n = 1 x 10 11 e/cm 2 is not caused by localization effects, 
impurity effects, etc., usually associated with the phe- 
nomenology of the 2-D metal-insulator transition. The 
enhancement of m*g* in the 1-valley 2DES of GaAs/AlAs 
systems is found to be strongly dependent on the spin 
polarizationQ, in strong contrast to the Si-MOSFET 
case. Our calculations [3| show that the physics of the 
1-valley system is dominated by the presence of a tran- 
sition to a fully polarized state which makes g* increase 
rapidly with r s as the transition density is approached. 
The 2-valley system shows no such transition to the spin- 
polarized state and is insensitive to the spin polarization. 

Given that perturbative many-body theory becomes 
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questionable for r s > 1, within the present context there 
are two approaches to evaluating the susceptibility en- 
hancement. The second derivative of the total free energy 
F(n, T) with respect to the spin-polarization £ gives a 
direct value for m*g* of the two-valley system. This re- 
quires the ^-dependent 2-valley energy which is not yet 
available from quantum Monte Carlo (QMC) simulations. 
However, we can evaluate F(n, £, T) using CHNC, and 
also show that the CHNC results agree with QMC data 
(available at ( = and T = 0). Another approach, which 
avoids the need for a full 4-component calculation is to 
build up the 2-valley susceptibility by noting that the 
one-valley energy F(n, C, T) is available at T = from 
QMC, and at any T ^ from CHNC. The coupling of 
the excitations of the two electron fluids in the two val- 
leys can be included in the coupled response function in 
a standard way. Then we find that the increase in g* 
in the GaAs 2DES is associated with the "blow-up" of 
the single-valley spin response, while the behaviour of 
the Si-MOSFET 2DES is related to the "blow-up" of the 
coupled-mode response. The static small-g limit of the 
spin response function provides the needed Xs/Xp- 

The objective of this paper is: (i) Present F(n, Q data 
for the two valley system by a full CHNC calculation 
involving the 10 pair distributions that are needed in 
the 2-valley system and establish the close agreement 
of the CHNC results with the QMC calculations, (ii) 
Construct the coupled-mode response functions using the 
well-established one-valley data and show that these re- 
sults are also validated by QMC and full CHNC results. 
That is, we present the CHNC results and compare with 
QMC data where available, and establish close agree- 
ment. We will not present detailed finite-T calculations 
(and hence m* calculations, as detailed in ref. @) in this 
paper, partly to limit the phase space that has to be 
studied, and since finite-T data are not presently avail- 
able from QMC for comparison. 

In section Mil we discuss the construction of the 
coupled-mode response functions which use only the one- 
valley exchange- correlation data to obtain the 2-valley 
behaviour via a physically motivated approximation to 
the inter-valley correlations. The compressibility pre- 
dicted via the small-g limit of the so constructed coupled- 
mode response is found to agree very well with that 
from QMC or the full 4-component CHNC energy cal- 
culation. This validates the couplcd-mode model used 
in the calculation (ref.|3j) of the m*g* enhancement in 
Si-MOSFETS. The full 2-valley energy calculations en- 
able us to examine the usual 1-valley local-density ap- 
proximation (LDA) in Si-MOSFETS and the corrections 
arising from coupled-mode effects. Finally we discuss the 
spin susceptibility enhancement obtained from these cal- 
culations, and the question of relating the electron-disk 
radius r s used in these calculations to the experimental 
densities. 



II. CHNC CALCULATIONS FOR THE 
2-VALLEY ELECTRON FLUID. 

We consider a 2-DES in a Si-MOSFET at a total den- 
sity n, with r 2 = l/itn, while the density in each valley 
v = a or b, is taken to be n v = n/2. Hence the r s param- 
eter in each valley becomes r sv = r s ^/2. Thus we do not 
consider density polrizations leading to n a ^ n\,. Also the 
electrons in both valleys have the same spin polarization 
£ and the same temperature T . This is consistent with 
recent studies which show that the valley splitting is very 
slight 0- If the two spin species be denoted by i = 1,2, 
we have a four-component 2-DES with 10 independent 
pair distribution functions (PDFs), viz., gij tVW (r). We 
define k = 1, 2 for the two spins in valley a, and k = 3, 4 
for the two spins in valley b, and write the PDFs as gu (r). 
The CHNC method for 2-DES has been described fully 
in ref. where the quantum fluid at T = is consid- 
ered to be equivalent to a classical fluid at a quantum 
temperature T q (r s ). Here, for the two valley system we 
use the same T q (r s ) as before, and embodies the essential 
"many-body" input to the problem. In brief outline, in 
CHNC we assume that the 2-D electrons are mapped on 
to a classical system where the distribution functions are 
given by a finite-T classical density functional form: 

9kl(r) = e~P {p ( r)Si i Svw+Vc ™( r)+Vc ( r '- [3kim (1) 

Here (3P(r) is a "Pauli exclusion potential" which acts 
only for parallel spins, i.e, if k — I. It is constructed such 
that gki{r) becomes identical with the non- interacting 
PDF, viz., yiAr) which is known from quantum mechan- 
ics when the Coulomb interaction V cou (r) and the associ- 
ated correlation corrections V c (r) are zero. The Coulomb 
interaction between two electrons in the equivalent clas- 
sical picture involves a correction arising from their mu- 
tual diffraction effects. Thus V cou (r) is obtained by solv- 
ing a two-electron Schrodinger equation. The result is 
parametrized by the form0, 

V cou (r) = (l/r)[l-«Bp(-fct h r)] (2) 
kth/k d Br = l.l587T? f 103 (3) 
k dBr = (2 7 rm*T c/ ) 1 / 2 , T 2 f = (T 2 +T 2 ) (4) 

Here kdBr is the de Broglie momentum of the scattering 
pair with the effective pair mass m* = 1/2, and T c f is the 
classical fluid temperature which reduces to T q at T = 0. 
The correlation potential V c (r) occurring in Eq.^is taken 
to be the sum of hyper-netted-chain diagrams inclusive of 
a bridge term. Thus V c is nonlocal and is a function of the 
gu (t) which have to be self-consistently calculated. The 
bridge term mimics the higher-order correlations which 
are not captured by the simplest HNC equations. These 
were shown to be important in 2-D electron systems in 
reference 0- Particles having identical indices (k = j) 
are restricted from close approach by the Pauli exclusion 
effect modeled by Pki(r). However, singlet-pairs of elec- 
trons, or electrons in two different valleys contribute to 
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strong Coulomb correlations, and hence a bridge term 
is included in all such "off-diagonal" PDFs. The bridge 
term Bki(r) for k ^ / applies to 6 different PDFs, and we 
have taken this to be given by the usual hard-disk func- 
tional form discussed in ref.[7j. (Totsuji and coworkers[8| 
have studied a more detailed implementation of the hard- 
disk bridge function in CHNC, while Bulutay et al.|9j 
have studied the 2-D CHNC without a bridge correc- 
tion). It should be emphasized that both the HNC ap- 
proximation as well as the need for a bridge function can 
be avoided by using the classical mapping to a quantum 
fluid (CMQF) where we use classical molecular dynam- 
ics (MD) to generate the PDFs etc., of the classical fluid 
under consideration. In such a scheme we use the pair- 
potential given by Eq.|2plus the Pauli potential in an MD 
simulation for a classical plasma at the temperature T c f. 
Such a CMQF-MD scheme would be numerically more 
demanding than the CHNC, much simpler than the full 
QMC simulations, and have the advantage of not making 
the HNC+bridge approximations. However, the 2-valley 
(4-component) system examined here has been studied 
by QMC and we use those results to confirm the validity 
of our methods. 

The main difference in the physics of the 1-valley sys- 
tem, and the 2-valley system arises from the preponder- 
ence of direct Coulomb interactions (from 6 PDFs in the 
2-valley, one in the 1-valley) over the exchange interac- 
tions (from 4 PDFs in the 2-valley, two PDFs in the 1- 
valley). This is the main reason for the lack of a tran- 
sition to a stable C = 1 state at low density. Since the 
transition to a ( = 1 state does not occur as r s increases, 
the g* remains insensitive to increasing r s , as found the- 
oretically (|3j) and experimentally [4|. 

The 10 coupled equations for gu(r) are self- 
consistently solved for many values of the coupling con- 
stant A applied to the Coulomb interaction, and the 
gkl(r : A) are used in the adiabatic connection formula 
to determine the exchange-correlation free energy of the 
2-valley 2DES. While our calculations are easily carried 
out for any value of £, T and r s , the 4-component QMC 
calculations at finite T, £ are a major computational un- 
dertaking which has not been attempted. However, at 
T = 0, C = 0, Conti and Senatore[l(j have presented 
QMC results for 2-D electron bilayers separated by a 
distance d^. They give total energies and also a fit to 
the correlation energy/electron, e c (r s ,C = 0,T = 0) at 
c?l = 0, i.e, the case where both electron gases reside in 
the same layer. In table H] we compare the 4-component 
CHNC with the available 4-component QMC data at 
d,L = 0. The energie s eJ r a ) QMC , are from the Rapisarda- 
Senatore nt-formula llll with the parameters quoted in 
Table I of Conti et alEj. 

These results show that the CHNC method provides 
a simple and accurate approach to the treatment of ex- 
change and correlation in the 4-component system. We 
exploit the simplicity of CHNC for calculations at finite- 
T and £ which are currently too expensive for QMC 
simulations. However, in situations where QMC re- 



TABLE I: Comparison of the total energy etot(r s ) and and 
the correlation energy e c (r s ), in atomic units at T = 0, £ = 
for the 2-valley 2DES obtained from CHNC, with the QMC 
data of Conti et al |lOI| . 

r a e tot (r s ) QMC e tot (r s ) CHNC e c (r a ) QMC e c (r 3 ) CHNC 



2.0 


-0.29302 


-0.29172 


-0.14315 


-0.14202 


10.0 


-0.08611 


-0.08647 


-0.04607 


-0.04649 


20.0 


-0.04641 


-0.04643 


-0.02577 


-0.02581 


30.0 


-0.03196 


-0.03183 


-0.01806 


-0.01795 



suits are available for the correlation energies, we adopt 
the QMC parametrizations. Thus one may use the 
parametrization of e c given by Conti and Senatore for 
the T — 0, C = case. For r s > 1 applications the 
Tanatar-Ceperlevfl^ form may also be used for the 2- 
valley system as well, with the parameter values ao = 
-0.40242, ai = 1.1319, a 2 = 1.3945, a 3 = 0.67883, fitted 
to a data base from r s =1 to 30. 

The T = 0, £ = 1 case is particularly interesting since 
this system (i.e, 2-valley system at density n) is math- 
ematically identical to the 1-valley system at the same 
density but with Q = 0, for the Hamiltonian considered 
by Conti et al., and by us in this study. In this case 
the 2-valley system has a two-fold symmetry since the 
energy is the same irrespective of the orientation of the 
spin in each valley. That is, £ = 1 means all the spins 
in valley a are oriented, while all the spins in valley b 
are also oriented, but independently of the orientation 
of the spins in a. This degeneracy would be resolved 
in real Si-MOSFETS but not in the model used here, 
or in Conti and Senatore. For instance, the three-body 
correlations for inter-valley interactions may be slightly 
different from those in the mira-valley interactions, and 
hence may require two different bridge parameters, to be 
determined variationally by an energy minimization us- 
ing the hard-disk reference fluid approach. We have not 
done this, and simply used the same bridge parameter as 
in ref. jjj for all interactions. In the QMC calculation this 
would require independent optimization of the model for 
back-flow corrections. Finally, the correlation energy of 
the fully spin-polarized (degenerate) 2-valley system can 
be parametrized using the Tanatar-Ceperley form with 
a = -0.19162, ai = 3.6123, a 2 = 1.9936, a 3 = 1.4714, 
in atomic units. 



A. The energy of unpolarized and polarized phases. 

The T — correlation energy at finite values of £ 
were calculated using the CHNC procedure and com- 
pared with the values predicted from the polarization 
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TABLE II: The correlation energy e c (r s , £) per electron, as a 
function of £, estimated using the 1-valley polarization factor 
of Eq. Q)J and from the full 2- valley CHNC calculation 





Jit 
c c 




CHNC 


Jit 
c c 


CHNC 




0.25 




0.25 


0.75 


0.75 


5.0 


-0.07686 




0.07757 


-0.06286 


-0.06434 


10.0 


-0.04518 




0.04562 


-0.03765 


-0.03831 


20.0 


-0.02529 




0.02535 


-0.02134 


-0.02151 


30.0 


-0.01772 




0.01764 


-0.01477 


-0.01504 


factor used for the 


one 


-valley 2DES. This has the form0 


P(r s ,0 


e c (r 8 ,0 


-ec(r„0) 


/-a(r s ) 
_ ^+ 


+ C (rs) - 2 


e c (r s ,l) 


-e c (r s ,0) 


2° 


<(r.) - 2 


a{r 8 ) 


= Cl- 


C 2 J 


'r s + C 3 /r 




l /3 (5) 



Here, £±=(1±£). It turns out that the coefficients C\— C4 
obtained for the one-valley 2DES, i.e., 1.5404, 0.030544, 
0.29621, and 0.23905 respectively, work quite well for the 
2- valley system as well, even at high r s . Thus, using the 
TC-type fit formulae for the 2-valley e(r s ,( = 1) and 
e c (r s ,C = 0), the estimated values of e c (r s ,() and the 
CHNC values are given in Table ILTI 

The finite-£ calculations show that there is no ferro- 
magnetic phase transition in this system at T = 0, since 
the total energy of the £ = phase is always more nega- 
tive than any polarized phase. This is expected from the 
dominance of the many off-diagonal terms contributing 
direct Coulomb interactions, but no exchange interac- 
tions. This was pointed out in ref. [j| where the insensi- 
tivity of m*g* to £ obtained from the theory was in excel- 
lent agreement with the experiments of Shashkin et al0. 
The stabilization energy AE of the £ = phase with re- 
spect to the fully polarized phase is 0.12734x 10~ 3 a.u. at 
r s =25, and diminishes to 0.89929xl0~ 4 a.u. at r s = 40. 
These are very small energy differences and within the 
error of the CHNC method and possibly of the 2-valley 
QMC calculations. However, the pattern of stability of 
the £ = phase holds for all r s investigated. Conti and 
Senatore|[l0j find a Wigner crystal phase at r s = 42 and 
claim that "in the intermediate regime the 4-component 
2-DES can be expected to mimic the 2DES ■ • • , with the 
appearance of a spin-polarized ground state" . We arrive 
at the opposite conclusion. 

B. LDA-type calculations for Si-MOSFETS. 

In most density-functional calculations of Si/SiC>2 
quantum wells, the Kohn-Sham exchange-correlation po- 
tential V xc (r) is calculated using the local density approx- 
imation (LDA) where the total density n(r) is considered 



TABLE III: Comparison of the exchange-correlation energy 
txc{r s ) per electron, and the Kohn-Sham potential V xc (r s ), in 
atomic units at T — 0, £ = for the 2-valley 2DES obtained 
from QMC-fit or CHNC, and from the LDA. 



r s 


t-xcij' s ) 


€-xc{t s ) ^ 


V xc (r s ) LDA 


VUrs) QMC 


2.0 


-0.38213 


-0.35535 


-0.55189 


-0.50271 


10.0 


-0.09007 


-0.08851 


-0.13133 


-0.12835 


20.0 


-0.04742 


-0.04699 


-0.06957 


-0.06874 


30.0 


-0.03241 


-0.03221 


-0.04769 


-0.04730 



without taking account of the valley degeneracy. In ef- 
fect, the electron gas is assumed to be a single electron 
gas at a density r s and its exchange-correlation energy 
£xc{rs) and the Kohn-Sham potential V xc (r s ) are calcu- 
lated at the given density. (We recall that V xc (r s ) is sim- 
ply the xc-contribution to the chemical potential /i xc )- 
Results for e xc and fi xc for electrons in a Si-MOSFET 
calculated correctly, i.e., taking account of its degener- 
ate valley structure, and in the usual LDA approach are 
compared in Table ITTT1 The full chemical potential fi, as 
well as the compressibility ratio K°/K calculated for the 
2-valley system, and that obtained within the LDA ap- 
proach, are also shown in Fig. ^ It is noteworthy that 
the total chemical potential has a minimum near r sm ~ 1. 
In effect, a low density electron fluid (r s > r sm ) whose 
chemical potential is equal to that of a high density gas 
(r s < r sm ) exists and this could lead to spontaneous den- 
sity inhomgenieties in these systems. 

When the actual electron densities in Si-MOSFETS are 
converted to effective r s units (see below), the r s range 
1-6 is the most important for device applications, and 
hence overestimates in V xc contained in the usual LDA 
approach could be significant. 



C. Pair-distribution functions in the Si-MOSFET 

system 

The PDFs, denoted by gki(r), embody the detailed 
particle correlations in the system. In Fig. [21 we display 
an illustrative set of pair-distribution functions. QMC- 
based PDFs have not been reported in the literature and 
hence we do not have a direct comparison. However, good 
agreement between CHNC and QMC- based PDFs have 
been found in other systems (e.g, 2DES, 3DES, and fluid 
hvdrogenl [7 HT3L HI ] . For £ = the four diagonal PDFs 
gkk are identical, and similarly, all the six off-diagonal 
PDFs are also identical. Hence there are actually only 
two distinct PDFs, just as in the single- valley case where 
gn and g\2 define the ( = case. These are shown for 
r s = 2, 5 and 10 in the top panel of Fig. [3 It is also clear 
that the correlation effects are mainly determined by the 
off-diagonal PDFs. At finite ( there are five independent 




FIG. 1: Left panels: Comparison of the total chemical po- 
tential (i i.e, the total Kohn-Sham potential, calculated at 
the total density n and r s , for the 2- valley system, and if 
the LDA were used (labeled 1-valley), ignoring the 2- valley 
nature. Right panels: Same for the compressibility ratio 

PDFs. There are two independent diagonal PDFs gn — 
<733 and 322 = 344- The three independent off-diagonals 
are 312 — 323 = 334 = 314, 313 and 1724- These are shown 
for the case r s = 10 and £ = 0.5 in the lower panel of 
Fig.H 



III. RESPONSE FUNCTIONS. 

In the following we discuss the linear response func- 
tions since the static small-fc limit can be related to the 
derivatives of the total energies that were calculated from 
CHNC or QMC if available. This enables us to verify 
a simple procedure for the construction of the 2-valley 
response functions using only single-valley exchange cor- 
relation dataQ- The density-density response function 
x(k, oj) will to be called the d-d response for brevity. It is 
expressed in terms of a reference "zeroth-order" Xnik, to) 
and a local-field factor (LFF), denoted by G(k,u). 

X (k, w) = X R(k, w)/[l -V k {l- G d (k, u)}x° R (k, u)} (6) 

The s-s response (or "spin susceptibility" ) is written as 

Xs (k,uj) = -n 2 B xl(k,u)/[l - V k {l - G s (k,u)}x° R (k,u)} 

(7) 



FIG. 2: Pair distribution functions of the two-valley system. 
For C — the 10 PDFs reduce to two indeependent PDFs. 
For C 7^ -0 there are 5 independent PDFs. 

where /is is the Bohr magneton. Note that our defini- 
tion of the spin-LFF differs somewhat from a commonly 
used definition 16]. Our form makes the d-d and s-s LFFs 
appear formally similar, at least at this stage. Hence a 
single discussion applies to both, and we drop the sub- 
scripts d and s. The reference function x R {k,w) has the 
form: 

Here k, q are two-dimensional vectors, while the corre- 
sponding single-particle energies are denoted by eg etc. 
The Fermi occupation number r may be chosen to 

be the non-interacting value, in which case x° is the 2- 
D Lindhard function 2] . An alternative choice is to use 
the fully-interacting density, evaluated from the fully- 
interacting chemical potential, as in density-functional 
theory (DFT). 

The small-A; limits of the local-field factors Gd(k,ui) 
and G s (k,ui) can be obtained from the second deriva- 
tives of the exchange-correlation free-energy functional 
Fxc{n, C) °f DFT, with respect to the charge densities 
or the spin densities^]- These second-order derivatives, 
together with the second derivative of F xc (n, £) with re- 
spect to T for a two-valley system were used in our m* 
and 3* calculations reported in ref.Q, using the CHNC 
technique. Here we look at the compressibility and 
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spin-susceptibility of the 2-valley system obtained from 
the small-fc limit of the coupled-mode response function 
(built up from 1-valley data) and compare it with that 
obtained directly from the 2-valley CHNC and QMC. 

A. Response functions of the two-valley system. 

The theory of the one-valley fluid can be used for the 
two-valley (4-component) fluid if there is no valley po- 
larization (i.e, the two valleys are assumed equivalent 
although distinct), as in ref.Q- As this may not be com- 
pletely clear from the abbreviated discussion in ref.Q, 
we present some details here. 

In the theory of classical fluids, the response functions 
are simply related to the structure factors, while the 
LFFs are simply related to the direct correlation func- 
tions of Ornstein-Zernikc(OZ) theory. Since this paper 
is directed more towards electron-fluid studies, we fol- 
low the language of the LFFS and the related response 
methodology rather than the OZ presentation. 

Let us indicate the species (which may be a valley index 
or a spin index) by u or v, taking the values 1,2, and let 
us consider a weak external potential (f> v (k, lu) which acts 
only on the electrons of species v. The external potential 
induces density deviations 5n v (k,uj) such that|l8|: 

5n u {k,u)) =y^Xuv(k,u)i<l>v(k ) uj) (9) 

V 

These equations define the linear d-d response functions 
involving the species u and v. The longitudinal dielectric 
function e(k,uj) is now given by 

l/e(k,u>) = l + V k J2xuv(k,u>) (10) 

u,v 

Here \% is the 2-D Coulomb interaction 2ir/k. Note that 
we are using effective atomic units (Hartrccs etc.) such 
that e 2 divided by the background dielectric constant is 
taken to be unity (see section lTV|l . To relate the response 
functions to the local fields, we consider the effective po- 
tentials U v (k,uj) such that: 

U v (k,u) = V k [l-G w (k,u)]5n v (k,w) (11) 

Thus the bare Coulomb interaction between the particles 
of type u and v is modified by the LFFs G uv . Hence 
we can write the density deviations 5n u (k,ui) in terms 
of the effective potentials and the zeroth order response 
functions as follows (we drop the k,uj labels for brevity). 

Sn u = xl[<t>u + V k J2( 1 - G uv)Sn v ] (12) 

V 

Now, by a comparison of the equations El an d 121 we can 
write down the response functions of the coupled two- 
component system in terms of the zeroth-order response 



functions and the LFFs. 







(13) 


Xl2 


= V kX lxl[l-G 1 2]/D 


(14) 


di 


= i_V fc x?[l-Gii] 


(15) 


di2 


= Vtaifl-Gia] 


(16) 


D 


= did 2 - di 2 d 2 i 


(17) 



We have suppressed the k,u> dependence in the above 
equations, and also not explicitly given X21, ^2 and efel 
for brevity. We now define a total coupled-mode response 
function XT(k, uS) via 

l/e(fc,w) = l + VfcXT(£,w) (18) 

Then the total two-component 2DES response is given 
by: 

Xt = [Xi+X2 + VkXixiG*}/D (19) 
Gs = Gn + G22 — G12 — G21 (20) 

In the following we some times denote xt by Xcm to em- 
phasize the coupled-mode nature of the total response. If 
we are dealing with a simple (one-valley) electron fluid, 
e.g., a partially spin-polarized electron gas, the species 
index v is simply the spin index. Notice that the cou- 
pling between the two systems (be they spins or valleys), 
replaces the individual denominators d± and c?2 by a new 
denominator D, common to both systems, and contain- 
ing the cross terms dij. That is, instead of the two sets of 
excitations given by the zeros of d\ and c?2, we now have 
a common set of "coupled-mode" excitations defined by 
the zeros of D. We emphasize that in this analysis we 
have not used any form of CHNC theory. 

All the response functions prior to the switching-on 
of the Coulomb interaction between the two valleys are 
known. The problem is to determine the cross terms dij, 
i.e, the inter- valley term, d uv , occurring in the coupled- 
mode denominator D, using only the free energy data for 
a single valley. If we consider the small-fc limit, we see 
that the terms dij occurring in D are directly related to 
the second density derivative or magnetization derivative 
of the free energy contributions Fij arising from the PDFs 
gij. We know these individual free energy contributions 
for the one- valley problem. 

Let us first consider the small-fc limit of the static re- 
sponse functions to make contact with the compressibil- 
ity and susceptibility sum rules. 

B. Small-fc limit of the static response. 

The small-fc behaviour of the static response is related 
to the second derivative with respect to the density (or 
magnetization) and this provides well known sum rules 
that we exploit here. For simplicity, and for comparison 
with the degenerate two-valley case, let us review the 
one- valley paramagnetic case £ = 0, and consider the 
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calculation Q of the small-/c, static (uj — 0) limit of the 
simple (one- valley) response function. 

Xv(n v ) = X°K)/[1 - V k (l - G vv ) X °(n v )}. (21) 

The density-density response function is associated with 
the proper polarization function n^. Dropping the 
species subscript v for the present, we have 



n = n°/(i + i4Gn ) 



(22) 
(23) 



n u = -x u 

The small-fc behaviour of this function states that 

n/n = k/k (24) 

The compressibility k is calculated via the chemical po- 
tential u, starting from the total free energy per unit 
volume obtained from the CHNC calculation. 

F = F + F x + F c (25) 
F = n v [el{n v ) + e v x {n v )] + ne c (n) (26) 



el = (l + C 2 )/2r s 2 
2^2 



(27) 



3^r s 
dF 

H = -r- = Mo + (J-x + Mc 
an 

1/K = n — — 

an 



[(l + C) 3/2 + (l-C) 3/2 ] (28) 

(29) 
(30) 



At T=0, the chemical potential is given (in Hartrees) by: 



n v ir - 



2(2/tt 



He 



(31) 



The compressibility calculated directly from the 4- 
component calculation should agree with that obtained 
from the coupled-valley formalism. There, the small-fc 
limit of the denominator of the density-density proper- 
polarization function is given by 

1 + V k G dX ° = k/ko 

Here Gd is the LFF of the density-density polarization 
function. Hence the denominator d±, or d2 of the d-d 
response occurring in EqEl for any particular species 
is available for the density-density response function in 
each valley. But the cross density-density LFFs, e.g., 
G\2, and the cross-denominators dn needed to form the 
coupled- valley forms are not yet specified. 

In the case of the spin susceptibility Xs , the role played 
by the compressibility is taken over by the spin-stiffness, 
which is the second derivative of the free energy F with 
respect to the spin polarization £. Here we have, for a 
single valley, 



Xs/xp = 1 + d 2 {r 2 sv F xc }/dC 2 



(32) 



Hence the denominators d u and d v of the spin- 
susceptibilities of the each 2DES are known, at the valley 



densities n v = n u = n/2, from a simple CHNC calcula- 
tion or from a QMC energy parametrization. However, 
here again the cross terms di2 and d%i, (which are equal) 
needed to complete the calculation of the coupled sus- 
ceptibility fEcis. n"3TT§|l are not yet specified. 

The cross term for the d-d response, or for the s-s 
response can be calculated if the free-energy contribu- 
tion F uv arising from the Coulomb interaction among the 
electrons in the two valleys is known. The interaction is 
among the electrons of valley u, at density n/2, and the 
electrons of valley v at density n/2. The inter- valley free- 
energy contribution F uv (n/2,n/2) is purely Coulombic, 
and hence it is clearly analogous to the correlation free- 
energy term arising from the antiparallel-spin PDF, i.e., 
9i2{n, C — 0) °f the simple one- valley 2DES at density n 
with two spin species. The case £ = ensures that the 
total density is split as n u = n v = n/2. Thus the d\2 
term needed for the spin-susceptibility calculation and 
the density-density response calculations are: 



s 
d 



■s d 12 = d 2 {r 2 F c [g 12 }}/dC 2 
d dia = -(2/n)d 2 F c [ gi 2}/dn 2 



(33) 
(34) 



Note that d\ is calculated from F xc (n/2,( = 0), while 
du is from F c (n,( = 0) of the simple one-valley 2DES. 
Hence, knowing d\, d\2, (which are equal to d%, and d%i 
since the valleys are degenerate), we can calculate the 
susceptibility enhancement Xs/xpi as we H the compress- 
ibility ratio k/kq of the interacting 2-valley 2DES, with- 
out actually solving the coupled system of 10 distribution 
functions needed in the full CHNC calculation of the 4- 
component system. 



Results for the compressibility ratio and the 
susceptibility ratio 



we consider the compressibility ratio Kq/K obtained 
by the coupled-mode analysis and from the 2-valley QMC 
data^Ji or equivalently, from the 4-component CHNC 
data, in Fig. The excellent agreement shows that the 
coupled-mode procedure for using the 1-valley data to 
generate 2-valley data is successful. 

The calculation of the susceptibility enhancement was 
given in ref. 3J. We consider the susceptibility enhance- 
ment in more detail. Equation 1321 involves the second 
£- derivative of the correlation energy. It is well known 
that the £- dependent QMC calculations are more prone 
to errors since a whole shell of spins need to be reversed. 
In fact, the 1-valley 2DES calculations of Rapisarda and 
Senatore^lJ using Diffusion Monte Carlo (DMC) predict 
a value of r s ~ 20 for the spin transition, different from 
that (r s ~ 26) predicted by Attacalite et al. 0] who also 
use very similar DMC methods. This is an indication of 
the type of uncertainty that may be had in ^-dependent 
QMC calculations. No C-dependent QMC data are avail- 
able for the 2-valley system. Using Eq. [3] we can write 
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FIG. 3: (a) Comparison of the compressibility ratio 
no/II cm = K I K obtained from the coupled-mode approach, 
and from the second-density derivative of the total free of the 
4-component system, (b) The energy estimate e c (rs,( — 0) 
from the coupled-mode form, and from the 4-component 
QMC0. 



an approximate explicit form at T = 0: 

e c (r s ,C) = e c (r s ,0) + (e c (r S! l)-e c (r S) 0))p(r s ,C) 
rf 2 e e (r s ,Q ABVl n , d 2 p(r s ,() 

de = As(1 ' 0) ^~ 

Thus the energy difference between the polarized phase 
and the unpolarized phase appears directly. This be- 
comes zero in systems like the one-valley 2-DES that 
show a spin transition. Even in the 2-valley system where 
there is no transition to a stable £ = 1 state, we can ex- 
pect the calculated Xs/xp or Xcm/xp to be quite sensi- 
tive to A_E(I, 0), and to the details of the form of the in- 
dependent function. We find that this is very much the 
case. In Fig.0]we display the coupled- mode evaluation of 
TO *3* = Xcm/xp with the value of Xs/xp obtained from 
the (^-second derivative of the energy obtained from the 
full 4-component calculation. We give curves labeled "4- 
component (a), and (b)", where the contribution from 
the d 2 e c /dC, 2 - derivative differs by ~5%. Clearly, this 
small change has a drastic effect on the m*g* evaluation. 
In the lower panel of Fig. [3] we compare the 2-valley cor- 
relation energy e c (r s , £ = 0) from the coupled-mode anal- 
ysis and from the direct 4-component QMC calculation. 
As expected, a small difference appears at low densities, 



FIG. 4: Comparison of the suscepribility enhancement 
Xcmjxp = m*g* obtained from the coupled-mode approach 
and from the second- £ derivative of the total free of the 
4-component system. The curves marked "4-components 
(a), (b)" differ by ~5% in the value of the second £- deriva- 
tive of the correlation energy, showing the strong sensitivity 
to this energy derivative. 

while the accuracy is good near r s ~ 5 where the rapid 
increase in m*g* occurs. This type of error is quite within 
the errors that are possible in the 4-component QMC or 
the 4-component CHNC, and in fitting to the polariza- 
tion dependence of the numerically calculated correlation 
energies. Coupled-mode formation implies that the exci- 
tation spectrum of the system no longer shows the fea- 
tures of the individual valleys, and hence is consistent 
with the conclusions of ref.[a| where no evidence for in- 
tervalley scattering was seen. 



IV. RELATION BETWEEN DENSITY n AND 
THE r s PARAMETER IN Si-MOSFETS. 

In the 2DES of GaAs/AlAs- based structures, the di- 
electric constants of the two materials are nearly iden- 
tical. The lattice constants are also well matched and 
hence the calculation of the effective atomic units needed 
in converting the experimental density n to the effective 
electron-disk radius (r s ) can be carried out in an unam- 
biguous, transparent way. This is important since many- 
body theory is formulated within the language of r s and 
effective atomic units. 
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Unlike in GaAs, the situation for Si-MOSFETS is 
more complicated. The lattice mismatch between crys- 
talline Si and most crystalline varieties (e.g., crystobalite) 
of SiC>2 turn out to be 35-40%. The dielectric constants 
are also strongly mismatched, being ~ 11.5 and ~ 3.9 for 
Si and Si02- The large lattice mismatch ensures that 
there is no sharp Si/SiOi2 interface. The reconstruc- 
tion of the Si atomic layers between the crystalline- Si 
and the SiC>2 bulk-like region still contain tetrahedral- 
bonding networks, but with strongly modified bond an- 
gles, bond lengths etc., characteristic of the amorphiza- 
tion of the Si layers immediately adjacent to the ox- 
ide layer. Many decades of experimental and theoret- 
ical work has gone into sharpening our understanding 
of this interface. More recently, first principles density- 
functional calculations by Carrier et al. [20j , starting from 
tight-binding models j^J have presented a clearer picture 
of the atomic arrangements near the interface region. A 
series of similar studies by Pasquerello et al. [23 establish 
the geometry of the Si/SiO^/vacuum interface. Thus a 
reliable atomic model of the Si/Si02 interface obtained 
via geometry optimization of the total energy is now 
available po|. The essential point is that the Si/SiC>2 in- 
terface contains approximately 5 regions containing crys- 
talline Si (c-Si), amorphized Si (a-Si), suboxide lay- 
ers, amorphized silicon dioxide (a-SzC^), and crystalline 
SiO-2 ■ These are indicated schematically below: 

[001] ^)|c-Si|a-Si|suboxides|a-SW 2 |c-SW 2 | (35) 

The amorphous (or bond distorted) regions of a-Si do not 
have a "conduction band" and should be considered as 
the true insulator that separates the 2DES which resides 
at the interface between c-Si and a-Si. Let the loca- 
tion of this amorphization edge be at z = z a . This edge 
can be defined to within a few atomic planes within the 
first-principles theoretical models. If we are dealing with 
a thick electron gas, then envelope- function methods for 
describing the form factor may be reasonable. Otherwise 
a more detailed atomic description involving Bloch func- 
tions is needed. In any case, if the electron gas is very 
thin, its growth-direction density profile may be consid- 
ered to be 8(z — z a ). That is, crystalline Si and a-Si 
flank the two sides of the 2-D electron layer, with a-Si 
playing the role of the insulator. The valence bonds of 
the a-Si still form a quasi-random tetrahedral network, 
even though distorted, and hence the "background" di- 
electric constant of a-Si is essentially that of Si That is, 
the effective dielectric constant: 



0.5(e s 







(36) 



often used for the 2DES in the MOSFET positioned at 
an abrupt Si/SiC>2 interface should be reconsidered. The 
second formula of Ando et al. (see appendix of ref.Q) 
for the conversion n to r s , using a mean e of 7.7 is not 
recommended. Instead, the first formula of ref. 2], i.e., 



r s /a* = 1.751 



10 12 cnr 



1/2 



11.5 



0.19m 



(37) 



is clearly the one consistent with the first-principles 
atomic structure of the Si/Si02 interface referred to 
above^. If we look at the Si-MOSFET literature, we 
find that the formula which uses the average dielectric 
constant of 7.7, appropriate for the abrupt Si/Si02 in- 
terface has been used by a number of authors. These au- 
thors use a value of r s increased by a factor of ~ 1.49 com- 
pared to what we recommend. Thus Pudalov et al.|23|. 
and Okamoto et al.[24j have used the mean dielectric 
constant of 7.7 for their calculation of r s . However, both 
these studies use the r s parameter mainly as a plotting 
variable in the figures, and not for any many-body calcu- 
lations. Hence a choice of r s which differs from that used 
in our work by a factor of ^1.49 is immaterial. In the re- 
view article by Kravchenko and SarachikjlJ, values of r s 
are further modified by the experimentally obtained m* 
to discuss the interactions in Si MOSFETS. Thus their r* 
is not used in the sense used in standard many-body the- 
ory. Hawang and Das Sarma 25] have also examined Si- 
MOSFET resistivities, using an impurity-scattering cal- 
culation which requires defining the effective background 
dielectric constant. They point out that their results are 
qualitative. Their results would not be affected by the 
choice of either formula given byref.Q, ie, using e —7.7 
or 11.5. Altshuler and Maslov[2^ actually consider the 
implications of the suboxide layer and how this could play 
a significant role in the theory. However, they too point 
out that their effort is essentially to indicate a "mecha- 
nism" rather than a theory of the metal-insulator physics 
of Sz-MOSFETS. Hence, once again, the results are too 
qualitative to make any difference. Similarly, the results 
of other workers [27| also do not discriminate sufficiently 
to make the choice of e a significant issue. 

Another class of problems where the choice of the av- 
erage dielectric constant is an issue is in calculating elec- 
tric subband energies |2j. Although the eigenvalues of the 
Kohn-Sham equation are not to be considered as effective 
excitation energies, such an assumption is often made. 
The input dielectric constant (which decides the effective 
r s ) enters into the exchange-correlation functions as well 
as the Poisson potential used. Most calculations are in 
the high density regime, and it turns out that, given the 
uncertainties of the quantum-well potentials and other 
parameters, the results can be equally well explained by 
a range of values of the effective dielectric constant. 

This situation becomes quite different when it comes 
to quantitative calculations for low-density MOSFETs, 
e.g., in the regime n = 1 X 10 11 electrons/cm 2 . Our 
CHNC calculations for m* and g* presented in ref.y] 
clearly favour the first formula, Eq. &\ of Ando et al. 
as the correct formula. This is also the formula that 
is consistent with the Car-Parrinello optimized atomic 
structure of the Si/SiC>2 interface obtained from the cal- 
culations of Carrier et al.[2(j and also of Pasquerello et 
a l-E2)- We believe that the problem of the correct dielec- 
tric constant at the Si/Si02 interface has received little 
scrutiny within the 2-D electron community in the past 
because there was no analytic many-body theory capa- 
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ble of giving quantitative results for low-density electron 
systems. Also, it is interesting to note that if the second 
formula of Ando et al. were used instead of the first for- 
mula that we recommend, then the calculated total r s is 
close to the n/2 value (per valley, ~ 1.414r s ) of r s cal- 
culated by the first formula. This is consistent with the 
calculations that we advocate, at the Hartee-Fock level 
(i.e, at the single-electron level). This fact can also lead 
to some confusion in assessing the validity of numerical 
calculations. 

The CHNC programs for electron gas calculations 
mentioned here and in ref. may be accessed via the 
internet at the address given in Ref. . 

V. CONCLUSION 

We can derive the following conclusions from this 
study. The CHNC method applied to a four-component 
electron fluid gives results in very close agreement with 



available Diffusion Monte Carlo calculations, without the 
use of any adjustable parameters specific to the 2-valley 
problem. The coupled-mode approach to constructing 
the 2-valley properties from 1-valley data is also fully 
confirmed. The calculation of the spin-susceptibility en- 
hancement x s /xp from the second £- derivative of the 
spin dependent energy is found to be very sensitive to the 
energy difference between the polarized and unpolarized 
phases and to the form of the polarization dependence. 
Nevertheless, the coupled-mode form is very successful in 
capturing the required physics. Thus the QMC, the 2- 
vallcy CHNC, and the coupled- mode approach based on 
the 1-valley data provide three independent methods for 
the study of the strongly coupled 2DES in Si MOSFETS. 
The three method are found to be in excellent agreement. 
Finally, we note that the methods used in this paper can 
be used to study electron/hole bilayers which are sepa- 
rated by a physical distance oIl, the present work being 
for d L = 0. 
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